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Abstract 

We consider a model for random deposition of monomers on a line with extrinsic 
precursor states. As the adsorbate coverage increases, the system develops non- 
trivial correlations due to the diffusion mediated deposition mechanism. In a nu- 
meric simulation, we study various quantities describing the evolution of the island 
structure. We propose a simple, self-consistent theory which incorporates pair cor- 
relations. The results for the correlations, island density number, average island size 
and probabilities of island nucleation, growth and coagulation show good agreement 
with the simulation data. 
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1 Introduction 



Models of random sequential adsorption (RSA) have been used to describe the pro- 
cess by which particles are irreversibly deposited without overlap onto a surface. 
They are relevant to studying the adsorption of gas molecules or colloid particles 
in solution || onto solid surfaces, or of large molecules on biological membranes ||. 
One of the key assumptions, is that the particles bind strongly to the substrate so 
that desorption and surface mobility are negligible on the time-scale of the experi- 
ment. Although providing a very simplified picture, RSA models have the virtue of 
being exactly solvable in dimension d — 1, for monomer and for fc-mer particles [|J. 

Many additional features have been added to this model in order to make it more 
physical, or to generalize its range of applications, and the literature in the field is 
vast The inclusion of diffusional relaxation on the substrate [H, [7[], leading to 
equilibration, or of cooperative effects, such as multi-site exclusion ||, have been 
examined. The possibilities that the particles are reflected back to the fluid or 
desorbed M, PL TO, or that the surface comprises more than one chemical species 



lTfl have also been considered. 

Another important possibility, motivated by many experimental studies of gas- 
metal surfaces [0, 13fl , is that if the particles lose just enough of their kinetic 
energy they can become trapped (physisorbed) in a mobile, temporary precursor 
state, from which they can be adsorbed (chemisorbed) at a far site at a later time. 
The mechanism of precursor mediated chemisorption, first postulated by Taylor and 
Langmuir ||14|| , was initially formulated as a statistical model by Kisliuk It was 
later adapted by King and coauthors [|H], [L3|, [16], [L| to include other effects, such as 



temperature dependence, desorption, molecular dissociation and pair interactions 
between the adatoms. Variations on the Kisliuk model have also been studied by 
other authors, both numerically and analytically || O, [TJ| O Most 
analytical treatments, however, are essentially mean-field like, in that they largely 
ignore the correlations between the state of occupation of different sites which arise 
through the precursor mediated deposition process, i.e. they assume that the rate of 
deposition at a given site only depends on the state of occupation of that site, and 
not on the neighbouring sites. 

In this paper, we consider a lattice model which is both a simplified version of 
the Kisliuk model, and a slight complication of the RSA model. Contrary to most 
previous studies, we examine the case where the mobility of physisorbed atoms is 
only possible on the top of occupied chemisorption sites (extrinsic precursor states), 
and where there is neither reflection nor desorption back to the gas phase (later we 
shall comment on how scattering and desorption can in principle be included into 
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the equations of motion). Each deposition attempt will, therefore, result either in 
the direct occupation of an empty site, or in the diffusion on the top of an island 
of occupied sites until the particle finds an empty site at the edge of the island, 
where it is irreversibly deposited. For simplicity, we also consider that the deposited 
particles are monomers, or atoms, in which case the lattice will eventually become 
full. As a consequence of the (extrinsic) precursor diffusion, the edges of islands are 
preferential sites for chemisorption and the growth of the larger islands is favoured. 
This introduces non-trivial correlations between sites, especially at late-times when 
diffusion is the dominant deposition mechanism, and makes the model unsolvable 
even for the simplest case of monomer deposition on a line. A model similar to the 



one adopted here has been examined by Becker and Ben-Shaul in 2d Their 
analysis, however, only applies to the early kinetics, as they treat the islands as 
uncorrected objects, completely neglecting the contribution of island coalescence to 
the growth process. 

Regarding the motivation for our study, we recall a somewhat unrealistic feature 
of the RSA model, namely that once a deposition attempt fails, the particle's posi- 
tion is 'randomized' before another adsorption attempt is made. As the transport of 
particles to the surface is diffusive, one might expect that a failed deposition attempt 
would be followed by another nearby adsorption attempt f2lfl . Since the intrinsic 
precursor seems to play no role in this effect, the present model may account for 
it, at least partially. Furthermore, it has been previously suggested by Cassuto and 
King [16||, that a kinetic deposition model without intrinsic precursor states and 
with negligible desorption, could well also explain the experimental data from some 
gas-solid systems, such as hydrogen on tungsten. Naively, one may also think that 
such a model could be able describe the slow deposition of liquid droplets on a plane 
surface, the spread of epidemics from immunized to non-immunized populations, or 
the growth of a forest where seeds have to be transported to an open field to find 
suitable conditions to develop. 

Although island formation and structure kinetics have been studied (see e.g. 
H ^Q, most work on RSA and precursor models, either analytical or simulation, 
has focused on determining the dependence of the coverage fraction 9 and sticking 
probability on the exposure, i.e. on the time t, and its value at the jamming limit 
(which is trivial only for monomers). In the present model, we assume that the dif- 
fusion time scale is small enough compared to the time scale of deposition attempts, 
and so the coverage is proportional to time or, equivalently, the sticking probability 
is 1. We examine the case of a substrate of dimension d = 1, and concentrate on 
studying the quantities describing the evolving morphology of the occupied regions. 
Namely, the pair correlation functions, the total number of islands, and the prob- 
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ability distribution of the island sizes and its moments. We have measured these 
quantities in numerical simulations, and developed a minimal theory incorporating 
correlations, to calculate them approximately. The results show good agreement, 
while at the same time differing considerably from those of the RSA model for 
monomer deposition in Id where there are no spatial correlations at all. 

A brief summary of this paper is the following. In section 2 we describe the model 
and define the formalism. The results of the numerical simulation are discussed in 
section 3. The theory is presented in section 4, and the results are compared with 
the simulation data. Section 5 is dedicated to some conclusions. 



2 Model 

We consider a one dimensional lattice with N sites, with N large enough to neglect 
boundary effects, i.e. we take N — > oo in the calculations. At each site i we define 
a variable Si, such that 

Si = if the site is empty 
= 1 if the site is occupied . 

Every time step 8t = 1 a site is chosen randomly and a monomer deposition attempt 
is made. If the site is empty the particle is adsorbed and the site becomes occupied 
irreversibly. Otherwise, the particle diffuses on top of the occupied region until 
it reaches an empty site where it is adsorbed (this is appropriately modelled by a 
random walk with traps [23 1). Only then is the next deposition attempted. The 



process repeats until the lattice is full. 

We assume here that the particle diffusion is rapid enough to be over before the 
time of the next deposition attempt, independently of the size of the island where 
it takes place. This may sound as an unrealistic assumption, especially if the island 
is of the order of the system size. We note, however, that due to the randomness 
of the process it is unlikely that a diffusing particle would interact or compete with 
the next particle to be deposited. This argument fails, of course, in the limit when 
the lattice is almost full (intermediate coagulation may occur) or in the case when 
diffusion is very slow (a gas of net precursor particles develops). 

Since each time step results in a deposition, either direct or mediated by diffusion, 
the time dependence of the adsorbate coverage, 6 = (Si), is trivial in this case: 

9(t) = t/N , < t < N . (1) 

One can easily incorporate in the model the possibility of scattering of a particle 
back into the gas phase at the instant of collision with the surface. If the scatter- 
ing probability is (1 — a), whatever the state of occupation of the site, the rate of 
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adsorption per unit time must be multiplied by a. We may then set a = 1 through 
a redefinition of the time scale. To account for the possibility that scattering may 
depend on the state of occupation of the site, we set the probability of no chemisorp- 
tion to be (1 — Da) if a particle lands on an occupied site. In this case, the rate of 
variation of the number of occupied sites is (setting a = 1) 

NdO/dt= 1(1 - 0) + DO = 1-(1-D)6 , (2) 

with initial condition 6(0) = 0. The solution is 

m = t-expHl-ZWM (3) 

This reduces to (jl]) for D = 1, and to the RSA behaviour for D = and for t <C N. 
When < D < 1 the full coverage 6 = 1 is attained at t = N\og(l/D)/(l - D), 
larger than N, due to the larger scattering from the occupied regions. We may also 
account for desorption, simply by adding an extra negative term (linear in 6) to 
the rate equation (|2]), effectively decreasing the value of D. The coefficient of the 
term may incorporate the probability of desorption from the physisorbed or from the 
chemisorbed states, or both. Allowing for desorption from the chemisorbed states 
would not only affect the coverage, however, but will also imply extra terms in the 
rates of growth and coalescence (eq. flT7D). For convenience we shall keep D / 1 
in the future expressions, as it allows us to distinguish the diffusion from the direct 
deposition terms. 

Next we define several quantities which will be useful to examine the morphology 
of the system. We shall call: 

• N(t) the total number of islands of occupied sites 

• N(L) = N(0L0) the number of islands with L sites, or number of rows of L 
consecutive occupied sites with at least one empty site on the left and right 

• iV(LOO) = iV(OLOO) the number of islands of L sites with at least two empty 
sites on the right 

• N(LOL') = N(0L0L'0) the number of islands of L sites separated by an empty 
site from an island of V sites on the right 



• iV(OO) the number of pairs of empty sites 
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• iV(OOO) the number of trios of empty sites . 

The generalization of the notation for more complicated configurations is obvious. 
The densities n associated with the above numbers, are defined by their ratio to 
the system size, e.g. n{t) = N(t)/N. Since the boundary effects can be neglected, 
the number of occupied regions (islands) is equal to the number of empty regions. 
On average the system has left and right symmetry, so iV(LOO) = iV(OOL) and 
N(LOL') = N(L'OL). 

For a Id substrate, the above quantities are easily defined in terms of the corre- 
lation functions of the local variables. Although the system is not in equilibrium we 
expect it to be translationally invariant, on average, due to the randomness of the 
deposition attempts. Using this property we have, for example: 

n(t) = £<(1- 3i)S i+1 ) IN = 9 - (SiS i+1 ) (4) 

i 

n(L) = Yl(^-Si)S i+1 ...S i+L (l-S i+L+1 ))lN 

i 

= ((l-S )S 1 ...S L (l-S L+1 )) (5) 
n(LOO) = ^((l-S i )S i+ i...S i+L (l-S i+L+1 )(l-S i+L+2 ))lN 

i 

= ((l-S )S 1 ...S L (l-S L+1 )(l-S L+2 )) (6) 
"(00) = £((1 - Si)(l - S i+1 ))/N = 1-29 + <£'o£i) (7) 

i 

n(OOO) + n(l) = £((1-^X1-^+2))/^ = 1-29 + (3 3 2 ) (8) 

i 

A number of normalization or hierarchical sum rules follow, some of which are: 

oo 

(9) 

oo 

= N(LOL') + iV(L00) (10) 

L'=l 

oo 

N(L00) + iV(000) (11) 

L=l 

N(0) = (1-9)N = VA^(OLO) = N(t) + N(00) (12) 

(13) 

It is easy to show that expressions (|) and @ (with n{L) replaced by @) are 
equivalent, using fixed or periodic boundary condition and the property Sf = Si. 
The probability of finding an island of size L and its moments, the average island 
size and its mean square deviation, are therefore: 

P(L) = N(L)/N(t) (14) 



N(t) 






L=l 


N(L) 


oo 




L'=0 


iV(00) 


oo 




L=0 


1-9)N 


oo 

= £^(0£0) 




L=0 


9N 


oo 
L=l 
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(L) = £ LP(L) = 9/n(t) (15) 



L=l 



{L2 \ L ,\ L)2 = i - *7 (»(*) £ ^)) ■ (is) 

There are three basic mechanisms by which these numbers may change with the 
deposition process (where A denotes variation): 

• nucleation: AiV(OOO) = 1 



• growth: AiV(LOO) = -1 ; AN(L + 1) = 1 

• coagulation: AN(LOL') = -1 ; AN(L + 1 + L') = 1 . 

The rates of occurrence of each of these events per unit time are: 
M000) , s 

Rn = = ^(000) 

00 1 J- nr/9 
R g = E N MLOO) + JV(OOL)] 

L=l iV 

oo 

= 2(n(00) - n(000)) + D Ln(LOO) (17) 

u^tl 1+fl( t +I,)/ V ') 

L=l L'=l iV 

oo 

= n(t) - n(00) + n(000) + D9-DY, Ln(LOO) 

L=l 

oo 

= 1 — (1 — D)6 — 2n(00) + n(000) -D^] Ln(LOO) , 

L=l 



where we have used (^)-(|T^). The factor L/2 accounts for a particle landing and 
diffusing on top of an island of size L towards the right (or left). For simplicity, we 
have assumed, and we will assume throughout, that in modeling the diffusion process 
the random walk with traps can be replaced by a random choice between right and 
left. We shall return to this point in section^. As expected, the sum of the rates ( |17|) 
yields the total rate of adsorption per unit time: J2i = 1 — (1 — D)6 = NdO/dt. 
The probabilities of nucleation, growth and coagulation per deposition event, P n , 
P g and P c , are then defined by the ratios: 

P, = Ri/(Nd6/dt) . (18) 

It is then straightforward to write the exact equation for the rate of variation of 
the number of islands, given by the difference between the rates of nucleation and 
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coagulation: 

dN(t) 

T. — — ti n — 1X C 

at 

oo 

= I- (D + 1)9 -2n(t) + DY, Ln(LQO) , (19) 

L=l 

where we have used fll2|) . Setting D = in (^) one obtains the RSA result n(t) = 
9(1—9). Using ( pT5|) and (§), it is possible to derive from (|I9|) an exact equation for the 
average island size (L). It is also illustrative to look at the growth of a single, isolated 
island. Neglecting coagulation, the equation for the island size reads: dL/dt = 
(2 + DL)/N, which solution (with L(0) = 0) is L(t) = (2/D)(exp(Dt/N) - 1). This, 
however, only gives the expected growth at early times. When island coagulation 
becomes important the growth of L should be faster than exponential, and as t — > iV 
(for D = 1) L should become of order N. 

3 Simulation 

We performed numerical simulations of the model described in section |2|. For conve- 
nience we used free boundary conditions, although the choice of boundary conditions 
should be irrelevant. The results were averaged over an ensemble comprising a great 
number (N s ) of system samples with different random deposition histories. For large 
coverages, the distribution P(L) and its second moment proved particularly sensi- 
tive to finite sampling and size effects. This is easily understood, as the size of the 
larger island at late-times (which controls the dynamics at this stage) can fluctuate 
by as much as N/2. The smoothness of the curves for the probabilities P n , P g and 
P c also depends on the number of samples used. We used N = N s = 50000, and 
verified that the systematic finite sampling and size errors were almost eliminated 
as we increased iV and N s up to these values. For ease of analytical treatment (eq. 
17) and for computational efficiency, we modelled the precursor layer diffusion with 



a random choice between right and left rather than with a random walk with traps. 
Although that mechanism does not take into account the starting point of the dif- 
fusion process, we expect this effect to be irrelevant because on average all sites on 
an island are equally likely to be chosen for a deposition attempt. 

3.1 Results 

The measurements of n{t), 1/ (L), ((L 2 ) - (L) )/ (L 2 ) and P(L) (at six different 
coverages) are displayed in Figures 1 to 4 (the pair correlations are shown in section 
4.1). We have divided the second moment by (L 2 ) rather then (L) 2 since the latter 



diverges. 
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It is interesting to observe the extent to which the RSA kinetics of deposition 
is modified by the additional diffusion mechanism. For comparison, we included in 
Figures 1 to 3 the plots (broken lines) of: 0(1 — 9), (1 — 9) and 9/(1 + 9), and in 
Figure 5 the plot of 9 L ~ l (l — 9), which correspond to the same quantities in the 
RSA model. We see from Figures 1 to 3, that the number of islands is smaller, 
and the average island size and island size fluctuations are larger with diffusion 
mediated deposition. This results from the increase in the rates of island growth 
and coagulation relative to the rate of island nucleation. Comparing Figures 4 and 5 
for the island size distribution, we can see a close agreement at small coverages, when 
all islands are still small. At intermediate coverages, the distribution spreads out 
(larger islands) in the diffusion case. At large coverages the difference is even greater. 
The RSA curve becomes uniform (all island sizes are equally probable), while the 
diffusion curve shows that the majority of islands are still small in size, but there is 
a minority of very large islands. At the late stages of diffusion mediated growth, the 
number of large islands is small and their size can fluctuate enormously, therefore 
the smoothness of the P(L) curve depends strongly on the sampling number. 

We have also examined the mechanisms responsible for the structure kinetics. 
Figure 6 shows the evolution of the probabilities of nucleation, growth and coagula- 
tion per deposition event, P n , P g and P c (eq. (0)). For comparison we have included 
the corresponding RSA curves (broken lines): P u ,rsa = (1 — 9) 2 , P 9y RSA = 29(1 — 9) 
and P c ,rsa = 9 2 (cf. (0)"(0) w hh D = 0). The most obvious difference occurs 
as 9 — > 1. In the RSA case the probability of growth vanishes. In the diffusion 
case, however, the probability of growth tends to a finite value 1 — P c , showing 
that a considerable number of gaps with two or more empty sites still exist at large 
coverages. 

4 Theory 

Due to the diffusion on top of the occupied regions, the probability of adsorption at 
an empty site depends on the state of occupation of its neighbouring sites and even 
of far located sites, if the site is at the edge of a large island. These correlations 
develop in the system as the coverage increases with time and diffusion mediated 
deposition becomes more likely. As a result, the density numbers defined in section |2] 
obey an infinite set of hierarchical coupled equations which cannot be solved exactly 
(an example of an equation in the top of the hierarchy is given by (|19D). One must, 
therefore, look for approximate solutions by truncating the hierarchy with some 
closure scheme. The simplest approximation consists in incorporating only the pair 
correlations in the system, which are then determined self-consistently. We shall 
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see that, by carefully choosing the multi-site functions decoupling scheme, such a 
simple approach, despite neglecting multiple correlations, is capable of capturing 
many qualitative and quantitative features of the system's behaviour. 

Let us denote the pair correlations between sites with different separations, as: 

p = = (3$+!) 

q = = (S t S l+2 ) (20) 

= (SiS i+n ) , n > 1 . 

We expect the correlations to decay with distance, i.e.: 

P>q> q {3) > ••• > 2 . (21) 

The quantities depending on the local variables of two sites only, can be expressed 
exactly in terms of the pair correlations. From ({§)-(§D we have: 

n{t) = 6-p (22) 
n (00) = l-26 + p (23) 
7i(000)+n(l) = 1-29 + q, (24) 

and the average island size ( |I~5| ) is given by 

(L) = jL. , (25 ) 

To write the multi-site functions approximately, in terms of the pair correlations, 
we decouple the higher order correlators into a product chain of pair correlators, each 
associated with an adjacent bond, as follows: 

/CO c o \ (S1S2) (S2S3) ... (S n -lS n ) (S n S n+m ) 

\D2/ ■■■ \ Vl/ Wn/ 

„n-l„(m) 

= qu i » (n>2,m>l) . (26) 

The normalization factors in the denominator assure that the RSA result is recovered 
in the decoupling limit. Then, from (§)-(§D we have: 

n(L) ~ ^(6-p)\ (L>1) (27) 
n(000) ~ 1 - 39 + 2p + q - p 2 /9 . (28) 

n(L00) ~ ^-l(e-p)[e 2 -9(p + q)+p 2 ] , (L>1) (29) 

Note that the limits of (n(00) - n(L)) and (n(000) - n(L00)) when L -> 0, which 
are {p — 9 2 )/p and (p — 9 2 + 9(q — p)) / p, respectively, although non-zero have a small 
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(31) 



value. Using these expressions and we can then write the deposition rates 

O as: 

Rn = l-39 + 2p + q-p 2 /9 

R 9 = (D + 2)(9-p) + 2(p 2 /9-q)+D9^- (30) 

u — p 

R c = Dp + q- p 2 /6 - D9 V -^- , 

6 — p 

and the mean square deviation of the island sizes (|16|) as: 

(L 2 ) - (L) 2 = _p_ 

(l 2 ) e + p' 

For the theory to be self-consistent, we must check that the expressions for the 
densities, although approximated are properly normalized. In fact, summing (p7|), 
with L from 1 to oo, gives 9 — p (cf. (||) and (f22"|)). Also, the sum J2T=i Ln(L) gives 
9 (cf. (0)). Hence, the probability distribution for the island sizes (0), obtained 
by dividing ( pT|) by n(t), 

P(L) = (p/9) L - 1 (l-p/9) , (32) 

is normalized to 1. This is a geometric distribution, as in the RSA case, but with 9 
replaced by p/9. Adding ( p8|) and (|29|) , with L from 1 to oo, gives 1 — 29 + p (cf. 
(|23|) and (|TT1) ). Other more complicated density numbers turn out to be consistently 
normalized too, as their sums yield the correct density number within the pair 
approximation. This is the case for n(LOL'), which satisfies the sum rule (|ToD. 

The next step of our approach is to determine the pair correlations self-consistently. 
There is an infinite hierarchy of coupled equations for the pair correlations, even 
within the pair approximation. The first two are the equations for p and q. We will 
now derive the exact form of these two equations. 

The p equation follows immediately from equation ([19]) for n(t), using ( p2[) and 



(|2|). It is more instructive, however, to write it down by inspection. pN is the 
average number of pairs 11, of two neighbouring occupied sites, which increases by 
one with island growth and by two with coalescence. Hence, using (|TTD, we have: 



,dp 
dt 



N^r = 2R C + R g 



2(D + l)9-2p-Dj2 Ln(LOO) . (33) 



L=l 



Within the pair approximation, using (^) and ([!]), and setting D = 1, we obtain: 

d9 =36 - p - e —p- (34) 
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The exact equation for g, although more complicated can also be written down by 
inspection. qN is the average number of pairs 1 — 1, of two occupied sites separated 
by a site in any state. It can increase by 1 or 2 with island nucleation, growth 
or coalescence. A careful consideration of all possibilities leads to the following 
equation: 

A f = 2 ± U(L0O0L>)+2 £ 1 + D ^ L + L ' )/2 N(L0L') (35) 

0,1 LL'=l iV LL'=2 iV 

00 1 _|_ DT' /O 

+ 2 £ N [N(LOOL') + N(L'OOL)} 

L=1,L'=2 ™ 

°° f 1 1 + D/2 1 

+ E { If [^(LOOOO) + iV(OOOOL)] + 1 [JV(L001) + iV(100L)] \ 



L=l 



+ 



£ \ 1 + N L/2 t^ 000 + N (W° L )} + 1 + D{1 n +L)/2 [JV(LOl) + iV(10L)]| . 



L=2 



Employing the hierarchical relations (P)- (^2|) and some obvious generalizations, this 
equation then simplifies to 

da 00 
N-j- = 2- 2(1 -D)6- (2 + D)n(l)-2n{000)-DY,L[n{L000) + n(L010)] . (36) 

^ L=l 

Evaluating the sum of L[n(L000) + n(L010)] within the pair approximation, using 
(|3) and (g§ and putting D = 1, yields 



_ = 2(6l+p- ff )-p/6l + -^-. (37) 

One can check that n(L000) and n(L010) are properly normalized within the pair 
approximation (their sums give n(000) and n(l), respectively). The approach also 
yields T,T=ol n ( Lom ) + n(L010)} = 1 - 29 + q, which is the exact result (cf. (f£f)). 

It is straightforward to derive the behaviour of p and q for small and for large 
coverage. As 9 — > the system is RSA like, so p, q and q^ should behave as 9 2 . 
Equations @ and © reduce top' = 39 + 0(9 2 ) and q' = 29 + 0(9 2 ), where primes 
indicate derivatives with respect to 9. Hence 

p = 3/2 9 2 + 0{9 3 ) (38) 
q = e 2 + 0{9 3 ) , (0->O) . (39) 

The factor 3/2 (confirmed by the simulations, Figure 7) shows that even in this 
regime there is an increase in the correlations relative to the RSA case. It results from 
the diffusion mediated growth of single site islands: with RSA there 2 possibilities 
for growth, and with diffusion there is a third one; therefore, there are 3/2 as many 
double site islands as in RSA. Let us now consider the limit when 9 — *■ 1. At 
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0=1,® and (0) yield the equations: (2 - p'(l))(l - p'(l)) = p'(l) - and 
(2 — g'(l))(l — p'(l)) = q'(l) — q^' (1). The first equation and inequality (|21~1) imply 
that 1 < p'(l) < q'(l) < q^'(l) < ... < 2. Hence, the solution is 

p'(l) = 2- (2-g'(l)) 1/2 = 2-(2-g (3) '(l)) 1/3 

q\l) = 2- (2-g (3) '(l)) 2/3 . (40) 

The actual value of the derivatives depends on the higher derivative qW , and thus 
on the truncation scheme adopted. The simulation results (Figure 7) suggest that 
p> (I) = q'(l) = 2, which implies (via (g0|) that g( 3 )'(l) = 2, and that 



p = q = g( 3 ) = 26 -1 + 0(1 -ef = e 2 + o{\-ef , (0^1). (41) 

This is consistent with the truncation schemes: q = 9 2 or q^ = 9 2 . 
4.1 Results 

The p equation (|34"D involves q; the q equation (|37|) involves q^ 3 \ etc. The nature 
of the approximation depends on how we close the hierarchy. In the simplest, first 
approximation we neglect the correlations beyond the nearest neighbours, i.e. we 
set q = qi = 6 2 , solve numerically equation ([34] ) for p = Pi, with initial condition 
p(0) = 0, and substitute p\ and q\ in the quantities of interest. In the second 
approximation we neglect the correlations beyond the second neighbours, i.e. we set 
q^ 1 — q^ — ® 2 1 an d solve the system of equations (|34])-(|3~T|) for p = p 2 and q = q 2 , 
with initial conditions p(0) = g(0) = 0. 

Figure 7 shows the simulation data, p s and q s , and the predictions from the first 
theory, pi and q\ = 9 2 . The correlations from the simulation decay with distance as 
in (|21~1), and the differences p s — 8 2 , q s — 9 2 are small but non-zero, as they should 
be since they establish the difference of behaviour relative to the RSA model. The 
agreement between p x and p s is quite good. As expected, it gets slightly worse as 9 
approaches 1 and the correlations between sites further apart become more relevant, 
but the correct asymptotic behaviour is obtained. Consequently, there is also good 
agreement in the results for the density of islands n(t) (Figure 8) and the average 
island size (L) (Figure 9), which are (exact) functions of p only. The agreement 
between q\ and q s (Figure 7) is, of course, less satisfactory. Figure 10 compares 
the theoretical and simulation plots for the probabilities of nucleation, growth and 
coagulation, P n , P g and P c (eq. ([3(]) with D — 1; cf. Figure 6). There is good quan- 
titative agreement up to 9 = 0.5, and there is still some qualitative agreement for 
larger coverages. The theory fails, however, to give the correct asymptotic behaviour 
as 9 — > 1: although there is a region where P g ~ 1 — P c , the theory gives P g — > as 
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9 — > 1. From fl3"U| ) we can see that P g ~ (1 — g')/ (1 — p') — 1 close to = 1. Hence 
the limit P g = is a consequence of having = which follows (via ( [40"1 ) 

from the truncation q = 6 2 . 

Nearest-neighbour correlations are sufficient to probe the presence of island 
boundaries. Since the above quantities depend essentially on island counting, their 
predictions are fairly accurate. Except in the early stages of deposition, however, 
long-range correlations need to be accounted for to correctly describe the spectrum 
of island sizes. Hence, as a result of the pair decoupling approximation (|26|), the 
theory predicts a geometric distribution for the island size probability (eq. fl32|)), 
the same as in the RSA case but with 9 replaced by p/9. The plots of P(L) are 
analogous to the RSA ones (Figure 5), though since p/9 > 9 each curve appears to 
correspond to a slightly larger coverage. The second moment of P{L) (eq. (|31])) is, 
of course, also incorrect: its limit as 9 — > 1 is 1/2 as in RSA, rather than 1 as in the 
simulation. From the plots of P(L) and its second moment, we find that the theory 
breaks down for these quantities for coverages over 0.2. 

Figure 11 compares the simulation data with the predictions from the second 

(3) 

theory, p 2 , q 2 and q 2 . As before, p 2 fits p s quite well, and q 2 fits q s even better. 
Consequently, the agreement in the results for n(t) (Figure 8) and 1/ (L) (Figure 
9) is also good: there is no major difference between the two theories, apart from 
the fact that now the curves lay below the simulation plots. A similar difference 
between theories is found (Figure 10) for the probability of nucleation P n = n(000), 
an approximate function of both p and q (eq. (^)). A worse agreement with the 
simulation than in the first theory is obtained, however, for the probabilities of 
growth and coagulation. We have also tested other plausible choices for the closure 
scheme, as the ones employed above are not unique, but found the results were either 
largely unchanged or incorrect. 

Finally, we note that p 2 lies over p s , while p\ lies below p s . q 2 , on the other 
hand, lies between q s and q\. Hence, the p correlations are underestimated in the 
first theory and overestimated in the second theory. This seems to indicate that the 
approach cannot be systematically improved by higher order truncations in the pair 
correlation hierarchy. This fact comes as no surprise given the uncontrolled nature 
of a self-consistent approach. Moreover, the results are more likely to be affected by 
the pair decoupling approximation (^) then by the order of truncation in the pair 
correlation hierarchy. 
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5 Conclusions 



We have studied numerically and analytically a simple, but non-trivial model for 
the deposition of monomers on a line. The particles can diffuse on the extrinsic 
precursor layer until they reach the edge of the island, where they are irreversibly 
deposited. As time progresses, island nucleation becomes less frequent, while the 
larger islands grow rapidly and merge with other islands. During this process, 
the precursor particles migrate for larger and larger distances with ever increasing 
probability, establishing correlations between sites further and further apart. As a 
result, the system develops a structure characterized by strong correlations whose 
range grows to the system size at full coverage. 

In the simulation, we looked at the evolving structure pattern by measuring the 
island density number, the island size probability distribution and its first and sec- 
ond moments, and some of the pair correlations. We also looked at the interplay 
between the direct and the diffusion mediated mechanisms of deposition, by measur- 
ing the nucleation, growth and coalescence probabilities per unit time. As expected, 
the results (Figs. 1-7) differ considerably from the RSA model, especially at large 
coverages, due to the increasing correlations. 

To explain and interpret these measurements, we proposed a simple, self-consistent 
theory which incorporates correlations to a minimum extent, i.e. local pair corre- 
lations. We considered two levels of approximation, depending on the closure of 
the hierarchical equations for the pair correlations. Altogether, the lowest level of 
approximation, accounting for nearest-neighbour correlations only, gave the best fit 
to the simulation data. The predictions of the theory (Figs. 7-9) proved very accu- 
rate for the nearest-neighbour correlator p, island density number n(t) and average 
island size (L). A fairly good prediction was also obtained for the probabilities of 
nucleation, growth and coagulation (Fig. 10). While the nearest-neighbour correla- 
tions are sufficient to distinguish occupied from non-occupied regions, the full range 
of correlations is required to distinguish the sizes of those regions. An accurate 
determination of the island size distribution P(L) and its second moment is, there- 
fore, beyond the scope this theory, and would, in principle, require the use of more 
complicated methods. 

We end with a comment on the 2c? systems, which are of more interest to ex- 
perimentalists. In this case, however, the same quantities are not easily expressed 
in terms of the local lattice variables, and one is faced with basic difficulties in the 
development of a useful formalism (most of section |2| would not be applicable) and 
in the analytical treatment. To illustrate the problem, we note that instead of the 
nearest-neighbour correlations, some non-local operator would be required to probe 
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the domain boundaries. It is desirable and possible, nonetheless, to perform nu- 
merical simulations of the 2d model, which would also enable the study of richer 
phenomena, such as percolating clusters. 
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Figure 1: Density number of islands, nit), as a function of the coverage, 6. Simula- 
tion data and RSA model. 




Figure 2: Inverse of the average island size, 1/ (L), as a function of the coverage, 6. 
Simulation data and RSA model. 




Figure 3: Mean square deviation of the average island size, 1 — (L) / (L 2 ) } as a 
function of the coverage 6. Simulation data and RSA model. 
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Figure 4: Probability distribution for the island sizes, P(L), at six different cover- 
ages, 9 = 0, 1 66, 0.333, 0.500, 0.666, 0.990, 0.999. Simulation data. 



hJ 04 

0.035 
0.03 

0.025 
0.02 

0.015 
0.01 

0.005 



: B -Q- B 



RSA 



= 0.1 66 



= 0.333 



W - = 0.500 



O - = 0.666 



= 0.990 



= 0.999 



B -Q- £]_ g 



-B- S-B-Qb-B-B- 



«-B-Q-fl- H ^ H _ H ^ & 



B-£>H-B^-a-[ 







A A A A A A A A A .A A A A & -A, A A A A A A A A A A A A A A A A A A A A A A A A 
I I I I T-^-.X >--!.. 1 T T T . . X < ?--<4>-(D-<j.-.rh-.*-A. A 1 I I I I I I I I I I I I T 







10 15 20 25 30 35 40 



Figure 5: Probability distribution for the island sizes, P(L), at six different cover- 
ages, 9 = 0, 1 66, 0.333, 0.500, 0.666, 0.990, 0.999. RSA model. 




Figure 6: Probabilities of nucleation, growth and coagulation, P n} P g and P c , as a 
function of the coverage 6. Simulation data and RSA model. 




Figure 7: Nearest-neighbour and second nearest neighbour pair correlations, p and 
q, as a function of the coverage 6. First theory (q = 6 2 ) and simulation data. 




Figure 8: Density number of islands, nit), as a function of the coverage, 6. First 
theory, second theory and simulation data. 




Figure 9: Inverse of the average island size, 1/ (L), as a function of the coverage, 6. 
First theory, second theory and simulation data. 




Figure 10: Probabilities of nucleation, growth and coagulation, P n} P g and P c , as 
a function of the coverage 6. Simulation data and RSA model. P n : hrst theory, 
second theory and simulation data. P g and P c : hrst theory and simulation data. 




Figure 11: Nearest-neighbour and second nearest neighbour pair correlations, p and 
q, as a function of the coverage 0. Second theory (q^ = 9 2 ) and simulation data. 



